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How to build a compact brane 
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This work deals with braneworld models in a five dimensional curved geometry with a single 
extra dimension of infinite extent. The investigation introduces a new family of models, generated 
from a source scalar field that supports kinklike structures described through the presence of a real 
parameter, capable of controlling the thickness of the warp factor that describes the hve dimensional 
geometry. The mechanism shows how to get a brane that engenders a compact profile. 
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I. INTRODUCTION 

Localized structures play important role in nonlinear 
science. In high energy physics, they may acquire dif¬ 
ferent features and appear as, for instance, kinks, vor¬ 
tices and monopoles, in one, two and three spatial dimen¬ 
sions, respectively mm- In the case of a single spatial 
dimension, kinks can be described by scalar field that 
self-interact through polynomial or nonpolynomial po¬ 
tentials. In the case of polynomial potentials, one usually 
considers fourth, sixth or higher-order polynomials, and 
in the case of nonpolynomial potentials, one has the sine- 
Gordon and the sinh-Gordon potentials as representative 
models. 

In the two-dimensional spacetime, topological struc¬ 
tures of the kinklike type can be used to study a diver¬ 
sity of problems in high energy physics, as we can see, for 
instance, in Refs. I3H7]. Also, they are of current interest 
to source braneworld scenarios in (4,1) spacetime dimen¬ 
sions, in a curved geometry with a single extra dimension 
of infinite extent [S{l6l. 

In this work, we start dealing with real scalar field 
in (I, I) spacetime dimensions, searching for kinklike 
structures and focusing mainly on the possibility to con¬ 
struct families of models, to be used to describe thick 
braneworld scenarios in (4,1) spacetime dimensions. The 
work is motivated from the recent investigation on kinks 
and compactons, which offered the possibility to generate 
thick brane with hybrid profile m Here, however, we 
introduce a new mechanism, capable of inducing compact 
profile to the warp factor that describe the geometry of 
the brane. 

The investigation concerns the study of scalar fields 
as source fields in a five-dimensional warped geometry, 
focusing on the possibility to navigate from models de¬ 
scribed by trigonometric interactions to other models. 
In order to get to the desired scenario, both in flat and 
in curved spacetime, we make use of the Jacobi elliptic 
functions, which are controlled by a single real parame¬ 
ter that connects trigonometric and hyperbolic functions. 
We start the investigation with the well-known standard 
(jA theory with spontaneous symmetry breaking, which 
supports kinklike structure, and we make good use of 
the deformation procedure introduced in [18] to obtain 
the model described by elliptic functions, which we can 


solve analytically, using the prescription of Ref. [T8| . This 
route is of interest since it offers the possibility of obtain¬ 
ing analytical results, which we then use to investigate 
the braneworld scenario. 


II. GENERALITIES 

We start with the Lagrange density 

C=\d^4^d^ct>-V{et>), ( 1 ) 

where 4> = represents the scalar field, p, = 0,1, 

and is the potential, which is used to identify the 

model under investigation. Here we are working with di¬ 
mensionless field and coordinates, for simplicity. The 
metric tensor is rj^i, = diag(l,—1), and the energy- 
momentum tensor has the form T^^’' = d^<j)d'^(j) — r]^''C. 
The equation of motion for the scalar field is given by 
dV/d(f> = 0; for static configuration we get 

d^^dV 
dx'^ d(j) 

We suppose that the potential V{(j)) can be written as 
V{(j}) = ^wl, (3) 

where = dw/d(j), with w = If we write the po¬ 

tential in the form above, we get to the first-order equa¬ 
tion 


which solves the equation of motion. The energy then 
becomes E= \w {(j){oo)) — w {(j) {—oo)))\, and the model 
can be seen as the bosonic portion of a supersymmetric 
model, with w = w{(j)) known as the superpotential. 

The strategy we shall follow in the current work is to 
use the deformation procedure put forward in [TH|. It 
allows to introduce new models from known models: we 
choose the model 0 and a specific potential V {(p) which 
we know to solve, to be the original model. We then 
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FIG. 1: Plots of the potential (181, for A = 0 (black, dotted 
line) 0.6 (green, dashed line), 0.9 (blue, dot-dashed line) and 
0.99 (solid, red line) 


FIG. 2: Plots of the solution (221 for some values of A, as in 
Fig. 1. 


introduce another field, X) and the function f{x), such 
that the new potential U(x) is given by 


u{x) 


f2 


( 5 ) 


derivative of the solution itself {r]o{x) « dx/dx) and other 
bound states, depending of the specific form of the model 
under investigation. This first-order framework help us 
to hnd analytical solutions in the curved five-dimensional 
spacetime which we are interested in, to construct the 
braneworld scenario. 


where = df/dx- The new theory is governed by the 
Lagrange density 


^X, d^x) = ^d^xd'^X -Uix), ( 6 ) 

and we have, for static field x = x{x) 


1 = -^' 


( 7 ) 


where W(x) is the superpotential of the new model. It 
is such that 




fix)) 

fx 


( 8 ) 


In this case, the energy of the static field conhguration is 


III. FAMILY OF MODELS 

An important model that supports kinklike solutions 
is the well-known model, described by the potential 

n^) = ^(i-<^")'. (12) 

This potential has the Z 2 symmetry and has minima at 
= ±1. The kinklike solution is 

(p(x) = tanh(a:). (13) 

We use this model and the deformation function given 
by the Jacobi elliptic sine 


A=|W(x(oo))-W(x(-oo))|. (9) 


/a(x) = sn(x, A), (14) 


As it was shown in [TB], if cj){x) is static solution of 
the original, (/)-field model, then the solution of the new 
X-field model is given by 

X{x) = /"^((/>(a;)). (10) 


where A is a real parameter, A € [0,1]; it is the modulus 
of the elliptic functions. We first recall that 

cn2(x,A)-bsn2(x,A) = 1, (15) 

dn^(x,A)-b Asn2(x,A) = 1. (16) 


Another important issue is that the new model is de¬ 
scribed by the superpotential that appears from (|^, so 
the kinklike structures also solve the first-order differen¬ 
tial equation 0- As one knows, the solution ( [l0| ) of the 
model 0 is linearly stable, with the stability potential 
having the form 


U{x) = W^^ + W^W^^^. ( 11 ) 

Here, x = xi^) is the solutio n ([To| of the first-order 
Eq. 0, for given by Eq. ([^. This potential en¬ 
traps the zero mode, r]Q{x), which is proportional to the 


The elliptic functions are interesting since they lead to 
both trigonometric and hyperbolic functions; for A = 0 
we have the usual trigonometric functions, and for A = 1 
we get to the hyperbolic functions: in fact, sn(x, 0) = 
sin(x), cn(x, 0) = cos(x), and sn(x, 1) = tanh(x) and 
cn(x, 1) = dn(x, 1) = sech(x). The functions sn(x, A) 
and cn(x, A) have period AK, while dn(x, A) has a smaller 
period, 2K, with 


^(l-t2)(l-At2)' 



( 17 ) 
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FIG. 3: The energy density (241 of the kinklike solution (221, 
depicted for some values of A, as in Figd] 


Thus, if we use the function (141 we get the potential 

1 


U(x,A) = -cd^(x,A), 


where cd(x, A) = cn(x, A)/dn(x, A). Also, 

W = -—In 

yX I dn(x,A) 


which leads to 


FFx = cd(x. A) = 


cn(x, A) 
dn(x,A)■ 


(18) 


(19) 


( 20 ) 


In Fig. we depict the potential (181. We see that 
the minima separate from each other as A increases from 
A = 0 to A = 1. We see that for A = 0, the potential is of 
the sine-Gordon type, and in the limit A —^ 1, it becomes 
constant, U —>■ 1/2. This is the behavior we want to 
use to build the braneworld scenario below. From the 
first-order equation one finds 


X' = cd (x, A). 

The solution for the above equation is given by 
x(x, A) = sn“^ (tanh(a;). A), 


( 21 ) 


( 22 ) 


which is depicted in Fig. for some values of A. The 
solution (221 is obtained from the deformation procedure 
m, and can be directly checked to be solution of ( [^ ; 
see Eq. (101. 

To understand the behavior of the solution in the limit 
A —)■ 1, we calculate its energy. The energy density Tqq = 
p{x) has the form 


P{x) = ^x'^ + U{x) = cd^ix, A) 

or, explicitly in terms of x, 

1 — tanh^(a;) 


p{x) = 


1 — A tanh (x) 


(23) 


(24) 


FIG. 4: The potential ( |44[ ) of the source scalar field in the 
braneworld model, depicted for A = 0 (black, dotted line), 
A = 1/3 (green, dashed line), 2/3 (blue, dot-dashed line) and 
A = 1 (red, solid line) 
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FIG. 5: The scalar field solution ( |45[ ) of the braneworld model 
(441, depicted for A as in Fig. 


The total energy is written as 


E(A) 


lAWl^^ln 


fl + VA] 

[i-VaJ- 


(25) 


It starts at Eg = 2, in the limit A —?► 0, and diverges in 
the limit A —>■ I. The issue here is that in the limit A —>■ 1, 
the potential becomes constant. As a consequence, the 
solution becomes a straight line, with constant derivative, 
leading the energy density to a constant value too, and 
making the energy to diverge as A approaches the unit 
value, as we illustrate in Figs. EJIllandlll 

The divergence found above makes the limit A —>■ 1 sin¬ 
gular in (1,1) spacetime dimensions. However, we iden¬ 
tify that the problem is in the potential, which becomes 
constant in the limit A —>■ 1. Thus, if we go to higher 
dimensions and add the model into a given, curved ge¬ 
ometry, the potential in the limit A —>■ 1 may describe a 
cosmological constant, circumventing the problem found 
in the (1,1) spacetime. 
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A. Braneworld scenario 

With the above motivation in mind, we now turn at¬ 
tention to the case of a braneworld model described by 
a five dimensional warped geometry, with a single ex¬ 
tra dimension of infinite extent IM- We consider the 
Einstein-Hilbert action which reads 

S = I d^x^\(^^R + Cy (26) 

Here we are using AttG^ = 1 and the Lagrange density of 
the source field is given by 

da^) = - U{^). (27) 

The metric is 

dsf = gabdx°‘dx^ = — dy^ (28) 

where ds\ = •q^^dx^dx'' and y describes the extra spatial 
dimension. The Einstein equation has the general form 

Gab = 2Tab, (29) 

As usual, we take A = A{y) and (f) = (j){y), and so we get 


6A'2 = (j)'^-2U, (30a) 

3A"-h6A'2 = -(f>'‘^-2U. (30b) 


Here prime stands for derivative with respect to y. We 
use (301 to get 


A" = 

We then introduce the hrst-order equations 


(31) 


0' = (32a) 

A' = -^W, (32b) 

which solve the equations of motion if the potential is 

given by 

= (33) 

An interesting feature that appears from the above 

braneworld scenario is that the total energy vanishes. To 
see this, we note that the energy density, which is given 
by 


Exp[A(y)] 



FIG. 6: The warp factor, depicted for A as in Fig. 


My) 



FIG. 7: The energy density (461, depicted for A as in Fig. 


B. Metric fluctuations 

Another general feature concerns stability of the grav¬ 
ity sector of the braneworld model. To show this, we 
redefine the y—coordinate to dy^ = and take 

the metric to the form 


ds^ = {qab + hab) dx°'dx^. (36) 


We then deal with conformal related metrics, i.e., gab = 
e^^gab- Using transverse-traceless gauge one can write 
G = —^dcd'^hab, which leads to the result 


GiV = -Idad'^hab + S 


daAdbA - dadbA + 


+ ^AXb + 5a&(9c9M + c>cA9M)]. (37) 


p{y) = 




+ U{cl>) , 


(34) 


(33), 


can be written in the form, with the use of Eqs. (32) and 

= ( 35 ) 


so it is such that the total energy usually adds to zero, 
since the warp factor exp (2A) vanishes exponentially, 
asymptotically. 


Thus, the /r^—components of the linearized tensor are 
+ -A'h'ij^^, — 35^1/ (a" -I- A'^) . (38) 

Here prime means derivative with respect to z. 

We can write the linearized energy-momentum tensor 
as 2T^]) = —ig^iv (A" -I- A'^). We then use = 2tI^J 
to get 

- + 3A'h'^^ = 0, (39) 













5 



ow 



FIG. 8: The energy density, depicted for A = 0 (black, dotted 
line), A = 0.3 (green, dashed line), A = 0.6 (blue, dot-dashed 
line) and A = 0.9 (p/3, red, solid line). 


FIG. 10: The potential of the new model, depicted for A = 0 
(black, dotted line), A = 0.3 (green, dashed line), A = 0.6 
(blue, dot-dashed line) and A = 0.9 {U /20, red, solid line). 
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FIG. 9: The zero mode, depicted for A = 0 (black, dotted 
line), A = 0.3 (green, dashed line), A = 0.6 {'x /2, blue, dot- 
dashed line) and A = 0.9 (xV20, red, solid line). 


which, under the transformation H^,y = e 
turns out to be 

(s, + (^-d, + (40) 

This equation has the form Sijj = and shows that 
the Hermitian operator S is non-negative, as it should 
be to describe linearly stable and consistent gravity sec¬ 
tor. The stability potential is written in the form 

U{z) = \a" + \a'\ (41) 

It has the usual profile and may support the zero mode 
(the graviton) as a stable state, bounded to the brane. 


C. Models 


Inspired from the study on the flat spacetime, we con¬ 
sider the function 


W = -^In 

Va 


1 — •\/Asn((/), A) 
dn((/. A) 


( 42 ) 


such that 


cn(^_ 

dn{(j), A) 

The potential then becomes 


(43) 


TTtj.-, cn(<(),A)2 1 2/I-yAsn(</),A)) 

^ 5a‘“ [- dn(».A) ) ^ 

It is depicted in Fig. for some values of A. The field 
solution is 


(j){y) = sn ^(tanh(2//2), A), 


(45) 


as given by Eq. (101, obtained from the deformation pro¬ 
cedure. It is depicted in Fig. for some values of A. 

Note that in the limit A —>■ 0, we have U{4>) = 
^(II cos(2(/)) — 5), and for A —1, one finds C/(^) = 
g — I ln^(cosh(^) — sinh(^)). The case A = 0 leads to a 
particular case of the sine-Gordon model studied in mi, 
for b = 2/3 and c = 1/2. 

We have been unable to hnd an analytical expression 
for the warp factor, so we depict the numerical solution 
for some values of A in Fig. Also, the energy density 
is given by 


piy) = 


- p2A(y) 


cn((/), A)^ 
4dn((/), A)^ 


In 


2 f l — x/~Xsn{(p^X) 
dn((p^X) 


3A 


g-difc 


It is depicted in Fig. [7] for some values of A. 


(46) 


One sees from Fig. ^ that if one increases the value 
of A, the amplitude of the solution also increases, but 
it makes no important change to its derivative near the 
origin. The effect of this in the geometry is to slowly 
narrow the warp factor, as it appears in Fig. 

The above results teach us that if we want to make 
the warp factor more localized, we have to increase the 
derivative of the scalar field solution near the origin. 
However, from the first-order equation (32a) we see that 


this can be implemented with the increasing of the value 
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Exp[2A(y)] 



FIG. 11: 
Fig.[Tg 


The kinklike solution (491, depicted for A as in 


FIG. 12: The warp factor, depicted for A as in Fig. |10[ 


of which is given by ( [4^ . We then make the fol¬ 
lowing change: we rewrite the superpotential of the de¬ 
formed model to the form 


W = - 


i-a^a 


1 I 1 — '/Asn(<^, A) 


dn((f>, A) 


(47) 


such that 


= 


1 cn(0,A) 

1 — A dn(((), A) 


(48) 


This modification requires that we exclude the case A = 
1, so we have to take A € [0,1) from now on. The new 
potential becomes = l7((/>)/(l — A)^, where I7((f>) is 
given by ( [44| , and the new solution is 

(^(y) = sn"^ (tanh(y/2(l - A), A). (49) 


We study this kinklike solution in (1,1) spacetime dimen¬ 
sions to see how it behaves as A increases in the interval 
[0,1). We depict the energy density and the correspond¬ 
ing (non normalized) zero mode in Figs. and re¬ 
spectively. They both show very clearly the tendency to 
concentrate around the origin as A increases toward unit, 
so it may contribute to make the brane compact. 

We then turn attention to the branewor ld s cenario. 
The new potential U(r/>) is depicted in Fig. 


10 


and the 

solution ( |49[ ) in Fig. [Hj We see that the solution be¬ 
comes thinner and thinner, and the amplitude increases 
as A increases in the interval A S [0,1). The new warp 
factor is depicted in Fig. and in Fig. we plot the 
energy density p for some values of A very close to unit. 
It is interesting to note that in the limit A —>■ 1, the kink¬ 
like solution behaves as the singular tachyon kink intro¬ 
duced sometime ago; see, e.g., Ref. m and references 
therein. However, in the case of A « 1 the model de¬ 
scribes braneworld configuration which engenders a reg¬ 
ular behavior, which is very well localized around the 
origin; as A increases toward unit, the warp function be¬ 
comes more and more localized, vanishing exponentially 
very strongly, and smoothly controlling the energy den¬ 
sity such that the total energy vanishes, as we illustrate 
in Fig. O 



-0.2 0 0.2 


FIG. 13: The energy density p, depicted for A = 0.93 (black, 
dotted line), A = 0.96 (blue, dot-dashed line, p/2) and A = 
0.99 {p/20, red, solid line). 


The behavior of the warp factor in Fig. 12 shows that 
its width decreases significantly as A increases toward 
unit, making the extra dimension compact. The proce¬ 
dure presented above offers a new braneworld scenario, 
which has nothing to do with the mechanisms explored 
before in m and more recently in m- 

To see if the gravity sector is robust, we focus on the 
stability potential (41) for the model under investigation. 
Here we follow }20j to study its profile, to see that it 
behaves smoothly even for A very close to unit. We then 
depict it in Fig. [^for some values of A close to unit. We 
note that the potential narrows as A increases toward 
unit, but it keeps the zero mode bounded to the brane. 

As A increases toward unit, the warp factor of the brane 
becomes thinner and thinner, but it is well different from 
the thin brane profile originally introduced in [9j. We 
have calculated the Kretschmann scalar to get 


= 40A'^ -b 16A"2 -b i2A'^A". 


(50) 


We checked that it behaves smoothly for A in the inter¬ 
val [0,1). It increases for increasing A, and diverges at 
the origin y = 0 for A = I, as expected. We then con¬ 
clude that the brane has a regular behavior as A varies 
in the interval € [0,1), and tends to become compact as 
A approaches unit. 
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U(z)) 



FIG. 14: The stability potential (41 1 for the model (47l, de¬ 
picted for A as in Fig. 13 with the same factors 1, 1/2, and 
1/20, respectively. 


IV. ENDING COMMENTS 

In this work we studied braneworld models described 
by a single real scalar field in a warped geometry with a 
single extra spatial dimension of infinite extent. We used 
the deformation procedure suggested in [T5] to build a 
scalar field model which is described by Jacobi elliptic 
functions, controlled by a single real parameter that ac¬ 
quires values in the interval A G [0,1). The interesting 


result is the presence of a new braneworld scenario, which 
describes a regular zero energy brane configuration with 
the warp function more and more localized around its 
center, as A increases toward unit. 

The braneworld model is controlled by a single real 
parameter A G [0,1), and its construction has nothing 
to do with the investigations described before in m 
As we can see from the energy density and the stability 
potential depicted in Figs. [T3| and [l4l the brane behaves 
regularly for A very close to unit. The model that we pro¬ 
posed leads to a new braneworld scenario, well different 
from the thin or thick brane profile introduced sometime 
ago. 

The proposed model is described by a parameter that 
controls the thickness of the brane in a very nice way, 
and may have applications of current interest to parti¬ 
cle phenomenology. An interesting line of investigation 
would be to study the presence of fermions and other 
fields, to verify how they can be entrapped inside the 
brane, controlled by A. 
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